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Abstract 

Efforts to extract information on magnitude and density depen- 
dence of the nuclear symmetry energy are discussed. The utilized data 
include those on mass dependence of the excitation energies to the 
isobaric analog states of ground states, as well as data on the diffu- 
sion of isospin in heavy-ion reactions. Results following from different 
observables are compared. 

The interest in nuclear symmetry energy has grown on account of the rel- 
atively recent availability of exotic beams, allowing to study systems with an 
increased range of relative asymmetry, 77 = (N — Z) / A, for a given mass num- 
ber A. Investigations have progressed focusing alternatively on the structure 
of nuclei, processes taking place in reactions and on collective excitations. 
Constraints on the symmetry energy would improve extrapolations from 
finite nuclei to neutron stars. 

Understanding of structural features, associated with the nuclear symme- 
try energy, can be advanced following the symmetries of nuclear interactions: 
charge symmetry, which is the symmetry under neutron-proton exchange, 
and charge invariance, which is the symmetry under rotations in neutron- 
proton space. In the context of charge symmetry, isoscalar quantities can 
be introduced, that do not change under neutron-proton interchange. If an 
isoscalar quantity F is expanded in r], the expansion contains even powers 
only: 

F{r,) = Fo + F27]^ + F47]^ + ... . (1) 
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Because of the lack of a linear term and because of small r/- values in nuclei, 
r] < 1/4, isoscalar quantities depend weakly on asymmetry. Examples of 
isoscalar quantities include nuclear energy and net nuclear density p = Pn + 
Pp. Isovector quantities are those that change sign under neutron-proton 
interchange. Example of isovector quantity is the neutron-proton density 
difference, pnp = Pn — Pp- An isovector quantity, expanded in rj, contains 
odd powers only: 

G(r?) = Gir? + G3r/3 + ... . (2) 

Notably, an isovector quantity divided by the asymmetry, G/rj, or by an- 
other isovector quantity, becomes an isoscalar quantity and, in this, weakly 
dependent on rj. 

The considerations above apply to quantities with a continuous depen- 
dence on rj. Microscopic, shell and pairing, effects can introduce discontin- 
uous changes. However, the considerations with expansion should apply to 
quantities that are averaged over the microscopic effects. In the context of 
charge invariance, the isoscalar quantities are those that do not change un- 
der rotations in neutron-proton space and the isovector quantities are those 
that rotate in a covariant manner. One consequence of the charge invariance 
of nuclear interactions is the appearance of isobaric analog multiplets across 
the isobaric chains, at about the same energy. Coulomb interactions break 
either symmetry, but may be accounted for in terms of correction terms. 

Up to the second order in rj, the nuclear contribution to the nuclear 
energy may be represented as 

E{N, Z) = Eo{A) + {N - Zf . (3) 

In simple mass formulas, the symmetry coefficient is usually assumed to 
be constant. However, the symmetry considerations alone do not prevent 
this coefficient from being ^-dependent, which we shall retain for generality. 
In ([3]), we may note a similarity to the energy of a capacitor in electrostatics, 
with capacitance C and charge Q, 

^ = ^0 + ^, (4) 

where Eq is the energy of the capacitor without charge. We can recog- 
nize that the asymmetry N — Z corresponds to the charge Q and A/{2aa) 
corresponds to the capacitance C. The analog of the capacitor voltage. 
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is the asymmetric chemical potential 

equal to, as in electrostatics, to the asymmetry charge divided by the ca- 
pacitance for asymmetry. 

Consistently with charge symmetry, the energy per nucleon in uniform 
matter may be represented, for low rj, as 

^{pn,Pp) = ^{p) + S{p)'n'^ . (7) 

However, microscopic calculations, such as [Ij, indicate that the r.h.s. of ([7]) 
represents the l.h.s. rather accurately all the way up to |??| = 1, for a wide 
range of p. In consequence, the two functions of net density, ^(p) and S{p), 
are sufficient to describe accurately the energy in uniform matter at different 
combinations of pn and pp. There is an interest in the expansion of those 
functions around the normal density pQ. By definition, Eq/A minimizes at 
po, but S generally has a finite slope, typically quantified in terms of the 
constant L: 

S{p) = al + \P^ + .... (8) 
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Here, = S{po) represents the symmetry coefficient for a large system 
dominated by normal density. Because of the minimum of Eq/A, the con- 
tribution from symmetry energy in ([7]) tends to dominate the pressure in 
neutron stars, at densities of the order of normal. 

As has been discussed, the net density, p = Pn — Pp, is isoscalar and, thus, 
should be weakly dependent on r], for a given A. We commonly parameterize 
p as 

1 + exp (^) 

with R = vqA^^^. On the other hand, the density difference, pn — Pp, is 
isovector. However, the ratio {pn — Pp)/il is isocalar. For a nucleus, the 
asymmetry r] in the latter normalization is a global quantity and it turns 
out to be more convenient to normalize the density difference in terms of 
the intense asymmetric chemical potential, to yield the asymmetric density: 

Pa{r) = ^[pn{r)-Pp{r)] . (10) 

Pa 

In the limit of a large system dominated by po, the two normalizations 
yield the same result. The two densities, p and pa both weakly depend on 



T]. Further, as will become partially apparent, the density pa is universally 
related to p. Out of those two densities, the densities of neutrons and protons 
in different nuclei may be constructed: 

PnAr) = l[pir)±^Pa{r)]. (11) 

The asymmetric density, that represents a formfactor for the isovector 
difference (p„ — pp), is related, on one hand, to the generalized symmetry 
coefficient aa{A) and, on the other, to a local value of the symmetry en- 
ergy S{p). The first relation follows from the fact that the capacitance for 
asymmetry may be represented as the ratio of asymmetry (charge in the 
electrostatic analogy) to the asymmetric potential (voltage), yielding 

aa{A) Pa J Pa a-a J 

The integral over pa is thus proportional to the capacitance for asymmetry. 
The farther the asymmetric density sticks out from the nuclear volume, the 
greater the capacitance for asymmetry. Otherwise, in uniform matter, we 
find from the definition of the chemical potential 

Pa = [Pn - Pp) , (13) 

P 

which yields for the asymmetric density 



Due to the short-range of nuclear interactions, the result ()14p is further 
expected to be approximately valid in weakly nonuniform matter. 

The anticipated weak dependence of the two densities on asymmetry is 
next tested in Fig. [TJ The separate panels show p (left) and pa (right), for 
different Skyrme interactions (from top to bottom), in Hartree-Fock calcu- 
lations of half-infinite nuclear matter [2] , at different asymmetries (different 
lines). In half-infinite matter, the shell effects are suppressed. In addition, 
the Coulomb interactions are switched off, eliminating the need for any 
Coulomb corrections. It is evident in Fig. [1] that the two densities indeed 
change very little with the asymmetry. 

Left panels of Fig. [2] next compare the isoscalar and isovector densities 
for the different interactions. Within the matter, the densities are fairly 
close to each other. However, in the surface area differences emerge that 




Figure 1: Net (isoscalar, on the left) and isovector (right) densities at different 
asymmetries in half-infinite nuclear matter, as a function of position, for different 
Skyrme interactions, after Ref. 



are strongly correlated to the L- value associated with the symmetry energy. 
The higher the L-value, the farther out is the isovector density displaced 
relative to the isoscalar density. This can be understood in terms of the 
behavior of the symmetry energy at low densities, seen in the right panel of 
Fig. [21 The isovector density pa follows the expectation of Eq. (fT4l) down 
to about the density of p Po/4 a representative classical return point. 
Low values of S{p), at p < po, enhance pa, as evident in Fig. [21 Densities 
p < po/4: are dominated by tunneling and local relations are not expected 
to hold. 

In the context of surface differences, for a large system we can get for 
the capacitance for asymmetry: 

2C . = ^ /d3.,.W = ^ Jd^rpir) ^^Jd^r (p. - p)ir) 

A_ 

(15) 



The last approximate equality follows from the fact that the two densities 




Figure 2: The left set of panels shows different densities in half-infinite nuclear- 
matter, as a function of position, for different Skyrme interactions, after Ref. [2]- 
The solid curves there represent the net density at asymmetry r] = 0. The short- 
dashed curves represent the isovector density at r/ = 0. Finally, the long-dashed 
curves represent expectations for the isovector density based on the local value of 
the symmetry energy. The right panel shows the dependence of symmetry energy 
on density, for different Skyrme interactions. 



are substantially different only in the surface region. We see here that the 
capacitance emerges here as a sum of two capacitances, one associated with 
nuclear interior, proportional to A, and one associated with the surface, pro- 
portional to A^/^. The surface capacitance, in terms of af , is tightly corre- 
lated with the slope-parameter L. Low values of L are associated with large 
surface capacitance and low values of af . Even tighter, and better physically 
justified, is the correlation between the ratios a^/of and L/a^ pj. Given 
the degree to which which Eq. (jl4p is followed for different interactions, the 
latter correlation is expected to be robust. 

Determination of the asymmetry coefficient with its mass dependence, 
by fitting nuclear masses with an energy formula, is difficult because the 
symmetry-energy contribution is small and its details compete against de- 
tails of other contributions to the energy [3]. However, upon generalizing the 
symmetry-energy term in an energy formula, effects of the symmetry term 
may be studied on a nucleus by nucleus basis, in isolation from other con- 
tributions to the energy of a nucleus. Specifically, we can observe that the 
symmetry energy may be represented in terms of the isospin of the nucleus 



(TjT^) as 

Ea = aa{A) ^^LZ^ =\aa{A) ^ . (16) 

This representation makes it apparent that, in the present form, the sym- 
metry energy is an isoscalar under charge symmetry, but not under charge 
invariance. However, if we replace the square of the third component of 
isospin with the square of the net isospin, we wih arrive at an isoscalar 
under charge invariance, required for the energy under that symmetry, 

Ea = 4a,(A) — = Aaa{A) ^ ^ ' 

This result should apply to a lowest state with a given isospin in the nucleus. 
Such excited states represent isobaric analog states (IAS) of the ground 
states of nuclei with a higher asymmetry in the specific isobaric chain. In the 
ground state, the quantum number for the net isospin agrees in magnitude 
with the third component, T = \Tz\. 

With the generalization, the excitation energy to an IAS becomes 

E^in) - E,{T,) =^{T2(T2 + 1) - Ti(ri + 1)} ^^g^ 

+Emic{T2,Tz) — -Ernie (?2, T2) , 

where we also account for corrections due to microscopic effects and defor- 
mation. We employ the corrections by Koura et al. [5] and obtain generalized 
symmetry coefficients on a nucleus-by-nucleus basis from 

AI\E 

using data on IAS compiled in Ref. |4J. Inverse values of the coefficient are 
plotted as a function of A~^^^ in Fig. [3l The coefficient values drop at low 
A, down to ~ lOMeV, which corresponds to an increase in capacitance per 
nucleon compared to heavy nuclei. In heavy nuclei, the coefficient values 
rise up to ~ 23MeV. 

Figure [3] demonstrates that the mass dependence of the symmetry coef- 
ficient may be fitted, at larger A, with 

a-\A) = {air' + {air' A-'l\ (20) 



which follows from Eq. (jl5p . This, principally, allows for a model- 
independent determination of the coefficients and af . With the correla- 
tion between a^/af and L/a}^, further, the value of L may be estimated. 
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Figure 3: Inverse values of symmetry coefficients plotted plotted vs inverse values 
of the cube root of mass number. Circles represent values derived from excitation 
energies to the states that are isobaric analogs of ground states when utilizing 
shell corrections by Koura et al. [5] . Line represents a linear fit to the results from 
analog states. Squares represent results from a Thomas-Fermi calculation [3]. 

Figure [3] shows also results from the simple Thomas-Fermi theory [3j that has 
the benefit of producing some curvature effect. Progressing in such a fash- 
ion, one can estimate the nuclear values of interest: = (31.5 — 33.5) MeV, 
af = (9.5 - 12) MeV and L ~ 95 MeV [2]. 

For the results such as above to be robust, though, the results following 
from an analysis need to agree with those underlying the theory, when re- 
alistic description of nuclei is applied. This is next tested in Fig. HI which 
shows sample results of calculations [6] of the symmetry coefficients within 
the Skyrme-Hartree-Fock calculations of unrealistically large [8j and realis- 
tic [?j masses, together with expectations based on the results for half-infinite 
matter [2]. When the formula (j20p is fitted to the medium and heavy nuclei 
found in nature, typically the coefficients result which are close to those 
established for the interaction in the infinite and half- infinite calculations. 
There are, however, some systematic differences. Namely, the fitted val- 
ues tends to be a bit higher and of values a bit lower than in the direct 
calculations. With this, the fitted values should lead to somewhat excessive 
L-values. For a few of the Skyrme interactions [2], however, we find large 
differences between the fitted and underlying values, such as for the SkI4 
interaction in Fig. El The situation is further illustrated in Fig. [6] that shows 
the relative deviations between the fitted and expected linear dependencies 
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Figure 4: Inverse values of symraetry coefficients, for the SkMP and SK272 interac- 
tions, plotted vs inverse values of the cube root of mass number. Symbols represent 
results from spherical Skyrme-Hartree-Fock calculations using the codes by Rein- 
hard for realistic [7, and unrealistically large 8^ nuclei with Coulomb interactions 
switched off. The lines represent, respectively, linear fits to nuclei in the mass region 
20 < A < 240 and expectations from calculations of half-infinite matter ^ . 
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Figure 5: Inverse values of symmetry coefficients, for the SkI4 interaction, plotted 
vs inverse values of the cube root of mass number. Symbols represent results 
from spherical Skyrme-Hartree-Fock calculations using the codes by Reinhard for 
realistic [7| and unrealistically large [8] nuclei with Coulomb interactions switched 
off. The lines represent, respectively, a linear fit to nuclei in the mass region 20 < 
A < 240 and expectation from calculations of half- infinite matter [2] . 



of l/fla on A^^^^, over 20 < A < 250 mass region, for about 150 different 
Skyrme parameterizations from the literature [2, 6j. For tlie majority of 
those interactions the deviations are small, but for a small percentage of 
those interactions the deviations are large. The interactions with the large 
deviations tend to have objectively unphysical features such as be unstable 
in the long-wavelength limit or exhibit unphysically strong nonlocality in the 
symmetry energy. Long-wavelength instabilities are signalled by one of the 
1 = Landau coefficients being lower than —1. Correspondingly, in the left 
panel of Fig. [6] the deviations between the fit and expectations are plotted 
against the lowest of the £ = Landau coefficients. It is seen that, in- 
deed, large negative values of the lowest £ = coefficient are associated with 
large deviations. Contribution of interactions to the nonlocality of symme- 
try energy, for Skyrme parameterizations, may be quantified in terms of the 
so-called coefficient in the term within the Skyrme Hamiltonian density [9] 



Excessively large magnitudes of ^ represent senselessly long range of inter- 
nucleon interactions. In the right panel of Fig. [6l the deviations are plotted 
against the ^-values. Again here it is seen that large C-values are associated 
with large deviations. It is apparent that before reliable corrections for 
curvature may be established, and firmer conclusions may be reached on 

, af and L, associated with the limits of infinite and half-infinite matter, 
the Skyrme interactions need to be filtered to reject those which exhibit 
significant nonphysical features of one type or another. 

Within central nuclear reactions, the symmetry energy and its density 
dependence can be assessed by examining the transport of neutron-proton 
asymmetry, or isospin, across the reacting system [lOpilj, see Fig. [71 Chang- 
ing symmetry energy can assist in the transport of isospin across the reaction 
zone |141ll5j. just as electric field can assist in the transport of charge along 
a wire. To minimize effects of the loss of isospin to the outside of the system, 
rather than transport across, different combinations of nuclei may be stud- 
ied in an experiment, symmetric and asymmetric (161 fTO| 113] . such as the 
neutron-rich ^^^Sn + ^^^Sn, more proton-rich ^^^Sn + ^^^Sn and mixed ^^''Sn 
+ ^^^Sn. The variable, which emphasizes the effects of isospin transport and 
facilitates comparisons between experiment and theory, is the ratio |16j 



^n-rich ^p-rich 

In the absence of isospin transport, the ratio R should reach values of ±1 
in the projectile and target regions. In the case of complete isospin mixing. 
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Figure 6: Relative rms deviation of the mass-dependent symmetry coefficients, 
expected from calculations of semi-infinite matter, from coefficients from spherical 
calculations, for a variety of Skyrme interactions, plotted vs a minimal (, — Q Landau 
coefficient (left panel) and vs the (^-parameter for the Skyrme energy functional. 




Figure 7: Stronger variation of symmetry energy with position enhances transport 
of neutron-proton asymmetry across the reacting system, sketched on the left. On 
the right, the results of transport calculations with different paramcterizations of 
the symmetry energy are compared to data |12) on isospin equilibration, from Sn 
+ Sn collisions at 50 MeV/nucleon, after Ref. [T3] . 
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Figure 8: Constraints on the parameters and L of the symmetry energy, from 
different indicated sources, after Ref. [13]. 

the expected ratio is ~ across the system. The convenience of the ratio 
in comparisons is that the asymmetry may be replaced by any quantity 
expected to be proportional to the asymmetry, with expected similar results 
for R. 

On the right of Fig. [71 a comparison may be seen of the data [12] on 
the ratio R, from Sn + Sn collisions at 50 MeV/nucleon, to the results of 
calculations within the ImQMD model [E]. In the calculations, symmetry 
energy of the form 

S{p) = 12.3 MeV {p/pof^ + 17.6 MeV (p/po)^' (23) 

has been used [13j- The value of , in the comparisons, minimizes in the 
vicinity of 7^ ^ 0.7. Implications for the symmetry energy from compar- 
ing central-reaction data to transport-model simulations are next shown in 
Fig. [8l Also results from the IAS analysis are shown there, together with 
results from analyzing collective excitations. While there is some level of 
convergence for conclusions on the symmetry energy, reached from different 
directions, there is no consistency yet. As has been mentioned, curvature 
corrections for teh IAS analysis are likely to lower both the L and values. 

To sum up, efforts to narrow down features of nuclear symmetry energy 
are advanced from different directions. We have mostly discussed here an 
effort from the structure direction and a little that from the reaction side. 



Qualitative understanding of connections between inputs and outcomes can 
be helpful in establishing the validity of the employed procedures as well 
as circumventing technical problems. To bring consistency between results 
obtained though different methodologies, it is necessary to understand sys- 
tematic errors associated with the methodologies. 

This work was supported by the National Science Foundation under 
Grant No. PHY-0800026. 
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